Spin-fluctuations and the peak-dip-hump feature in the photoemission spectrum of actinides 
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We present first-principles multiband spin susceptibility calculations within the random-phase approxima- 
tion for four isostructural superconducting PuCoIns, PuCoGas, PuRhGas, and nonsuperconducting UCoGas 
actinides. The results show that a strong peak in the spin-fluctuation dressed self-energy is present around 0.5 
eV in all materials, which is mostly created by 5/ electrons. These fluctuations couple to the single-particle 
spectrum and give rise to a peak-dip-hump feature, characteristic of the coexistence of itinerant and localized 
electronic states. Results are in quantitative agreement with photoemission spectra. Finally, we show that the 
studied actinides can be understood within the rigid-band filling approach, in which the spin-fluctuation coupling 
constant follows the same materials dependence as the superconducting transition temperature Tc. 
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The discovery of superconductivity in PuCoGas 111 and 
soon thereafter in isostructural PuRhGas |2 1, and PuCoIns O 
(collectively called Pu-115 series) has revitalized the inter- 
est in the spin- fluctuation mechanism of high-temperature su- 
perconductivity. In particular, a systematic study of spin- 
fluctuation temperature versus superconducting transition 
temperature Tc indicates that Pu-115 compounds lie in be- 
tween the Ce-based 4/-electron heavy-fermion and (i-electron 
superconductors (cuprates and pnictides) |4 |. Within the ac- 
tinide series, the duality of correlation effects in plutonium 
compounds stems from Pu's position between the itinerant 5/ 
states of uranium | 5 | and the localized 5/ states of ameri- 
cium f6^ . This makes Pu a unique candidate to define the in- 
termediate coupling regime of Coulomb interaction in which 
neither the purely itinerant mean-field theory nor the strong- 
coupling Kondo lattice model hold exactly — a prototypical 
example of strongly correlated electron systems |7|. On the 
other hand, the diagrammatic perturbation theory of fluctua- 
tions can still be applied as long as the Hubbard U ^ W, 
where W is the non-interacting bandwidth 1 8 1. Therefore, it is 
important to characterize the evolution of the spin- fluctuation 
excitations in Pu-1 15 compounds, which will help to delineate 
the role of spin-fluctuation mediated superconductivity in /- 
electron systems. 

Photoemission spectroscopy (PES) has revealed a strong 
spectral weight redistribution in the single-particle spectrum 
with a prominent peak-dip-hump structure around 0.5 eV in 
PuCoGas |9|. This feature has been interpreted as the sep- 
aration between itinerant (peak) and localized (hump) elec- 
tronic states of the 5/ electrons dOlIID- To provide insights 
into this PES structure, we present a first-principles multi- 
band spin susceptibility calculation within the random-phase 
approximation (RPA). The results show that a considerably 
large amount of the spin-fluctuation instability is present in 
the 0.5 eV energy range which originates from the particle- 
hole channel between 5/ states. The resulting self-energy cor- 
rection due to spin fluctuations is calculated within the GW 
approach, which quantitatively reproduces the observed peak- 
dip-hump PES feature in PuCoGas . 

We interpret the spin- fluctuation effects on PES along the 



same line as the localized-itinerant duality discussed above. 
The fluctuation spectrum creates a dip in the single-particle 
excitations due to strong scattering. The lost spectral weight 
(dip) is distributed partially to the renormalized itinerant states 
at the Fermi level (peak), as well as to the strongly local- 
ized incoherent states at higher energy (hump). The coher- 
ent states at the Fermi level can still be characterized as Bloch 
waves, though renormalized, whereas the incoherent electrons 
are localized in real space exhibiting the dispersionless hump 
feature. We perform calculations for the actinide materials 
PuCoIns (Tc = 2.4 K), PuCoGag (Tc = 18.5 K), PuRhGas 
(Tc = 8.7 K), and UCoGas (Tc = K), which show that 
the computed spin fluctuations play a significant role for the 
systematic evolution of the electronic band renormalization 
and spectral weight redistribution across these compounds. 
We also deduce the computed spin-fluctuation coupling con- 
stant A, which follows Tc as we move across the series from 
PuCoIns^PuCoGas^PuRhGas^UCoGas, suggesting that 
spin fluctuations play a crucial role in the pairing mechanism. 
The results also demonstrate that the actinides can be under- 
stood within a unified description of rigid-band shift of the 5/ 
electrons close to the Fermi level (hole doping). 

Intermediate coupling model— We calculate materials spe- 
cific first-principles electronic band structures, including spin- 
orbit coupling, within the framework of density functional 
theory in the generalized-gradient approximation (GGA) 1 12 |. 
We use the full-potential linearized augmented plane wave 
method of WIEN2k 1 13 1. The calculation is performed with 
40 bands to capture the ±10 eV energy window of relevance 
around the Fermi level. The non-interacting susceptibility 
in the particle-hole channel represents joint density of states 
(JDOS), which can be calculated by convoluting the multior- 
bital Green's function Gsp(k, zcj^) (s^p are orbital indices), 
to obtain (spin and charge bare susceptibility are the same in 
the paramagnetic ground state)| 14 |: 

X%qMi^) = -J^^Gsp(k,iuJn)Gqr(k^q,iuJn^^^) 
k,n 

Within the RPA, spin and charge channels become decoupled. 
(We ignore particle-particle as well as weaker charge fluctua- 



tion processes.) In the spin-channel, the collective many-body 
corrections of the spin fluctuation spectrum can be written 
in matrix representation: x = X^[l ~ ^sX^]~^- The inter- 
action matrix Us is defined in the same basis consisting of 
intra-orbital /7, inter-orbital V, Hund's coupling J and pair- 
scattering J' terms. |[T4lfT6ll In the present calculation, we ne- 
glect the orbital overlap of eigenstates, and hence becomes 
a diagonal matrix and J = = 0. 

Using the GW approximation, where G represents the 
Green's function and W is the interaction vertex, we write 
the spin-fluctuation interaction vertex following Ref. |17J 

as Vpqrsi^.^^) — \UsX" {q^^)Us ' The Feynmann- 

L J pqrs 

Dyson equation for the imaginary part of the self-energy in 
a multiband system with N sites at T = is (for details see 
the supplement [18 J ) 



s 

(^^)>^iV,,(c^-^^), (2) 



for > 0, where the density of states is given by A^rs(e) = 
— Im[G^5(/c, e)] /tt. (For cj < 0, the only changes are 
that the upper limit of the integral is and the argument of 
Nrs is ^ — which is < 0.) F is the vertex correction 
discussed later. For a more accurate calculation, one needs 
to account for the anisotropy in F(q, Vt). In the present case, 
where the spin-fluctuation spectrum is considerably isotropic 
(see Fig. [2]), it is justified to use a momentum-averaged spin- 
fluctuation function. 

We use Eq. (2) to compute the imaginary part of the self- 
energy from the first-principles band structure. The real part 
of the self-energy, Tl^^{uj\ is obtained by using the Kramers- 
Kronig relationship. Finally, the self-energy dressed quasi- 
particle spectrum is determined by Dyson's equation: = 



S. The full self-consistency in the GW approxima- 



tion requires the dressed Green's function G to be used in x^. 
This procedure is numerically expensive, especially for multi- 
band systems. To overcome this burden, we adopt a modi- 
fied self-consistency scheme, where we expand the real part of 
the self-energy E^^ ^ {1 — Z~^)uj in the low-energy region 
where T^^q 0. The resulting self-energy dressed Green's 
function is used in Eqs. (l)-(2) which keeps the formalism un- 
changed with respect to the renormalized band f ^ = . In 
this approximation the vertex correction in Eq. (2) simplifies 
to r = 1/Z according to the Ward identity. We note that all 
calculations are performed by solving matrix equations, while 
the results shown below are for the trace of each quantity. For 
brevity, we drop the symbol 'trace' altogether. 

Results— Figure [TJpresents the calculated GGA band struc- 
ture in (a) and corresponding non-interacting DOS in (b) for 
all four materials studied here. Notice that the low-energy 
band structure remains very much the same for all materi- 
als. It only shifts upward in energy in moving along the se- 
ries PuCoIn5^PuCoGa5^PuRhGa5^UCoGa5. This behav- 
ior can be accounted for by a rigid band shift, see insets to 
Fig. 1. The Pu-115 compounds show two sharp peaks in the 
DOS just below and above Ef, which are mainly originated 
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FIG. 1. (Color online) (a) First-principles GGA electronic band- 
structure calculations for various Pu-115 and UCoGas actinides near 
Ef- (b) Corresponding DOS in the low-energy region of present in- 
terest. The arrows mark the relevant particle-hole excitations. Insets: 
Low-energy regions of dispersion and DOS showing that all materi- 
als are related by a rigid shift of bands in this energy scale. 



from the 5/ electrons of Pu atoms. The 3d (or M) and Ap 
(or 5p) electrons of the reservoir elements are not important 
in this energy scale [See Ref s . [TTI and [T9l for partial DOS]. As 
the DOS at Ef decreases in going to UCoGas (see cyan lines 
in Fig.[T]), most of the 5/ states move above Ef, reducing the 
correlation strength to a large extent. 

Projections of the computed spin- fluctuation vertex, 
F(q, cj), are plotted in Fig.[2]as a function of excitation energy 
along the high- symmetry momentum directions. Our choice 
of the screened Coulomb term U satisfies the intermediate 
coupling approximation of U/W ^ 1. As seen from the band 
structures in Fig. 1(a), the average bare bandwidth for all ma- 
terials near the Fermi level is of order of 1 eV. Hence, we set 
U=l qV for all compounds, which is below the critical value 
of a magnetic instability, that is, Ux^{q,^ = 0) < 1 for all q. 
Note that our screened U for the spin-fluctuation calculation 
is smaller than that used in LDA+/7 type calculations, where 
a rather large value of U 3 eV was introduced into the local 
orbital basis ll2QH22ll . 

All spectra split mainly into two energy scales (at higher en- 



3 



PuCoIiij PuCoGa^ 



PuRhGa; UCoGag 



0.5- 
> 

a 






r Mxrz 



J7 log<V> 



(e) 



■ PuCoIn5 PuRhGa5 
' PuCoGa5 UCoGa5 




Excitation energy, Q [eV] 



FIG. 2. (Color online) The spin-fluctuation vertex F (q, ^) is plot- 
ted along high- symmetry directions in (a)-(d). Panel (e): The cor- 
responding (y{Q))^ averaged over 3D momentum space. All the 
calculations are performed for ±10 eV energy window, but results 
are shown only in the relevant energy region. 



ergy, no other prominent peak is seen in the computed spectra 
up to 10 eV and thus not shown). Corresponding momentum- 
averaged values (V)^ are fairly similar for all Pu-115 com- 
pounds, but notable different for UCoGas. The low-energy 
peak arises from the transition between the 5/ states just be- 
low to above Ep (within the RPA, the peak shifts to lower- 
energy), see gold arrow in Figs. 1(b) and 2(e). The high en- 
ergy hump comes mostly from the transition of the second 
peak in the DOS below (hybridized d- and p- states also 
contribute 1 19|) to the 5/ states above Ep sls marked by the 
magenta arrow in Figs. 1(b) and 2(e). For UCoGas most of 
the 5/ states shift above Ef and thus intra-orbital spin fluctu- 
ations do not survive, while the inter-orbital spin fluctuations 
move to higher energy. 

The coupling of the spin fluctuations to the quasiparticle ex- 
citations gives the self-energy correction in Eq. 2. The imag- 
inary and real part of S are plotted in Figs.jSja) and[3jb), re- 
spectively. Note that T>"{uj) shows a peak-dip-hump feature, 
although strongly enhanced by the DOS in comparison with 
(y{rt))^. Both the low- and high-energy features move to- 
ward = as the 5/ states shift toward Ef across the series 
PuCoIn5^PuCoGa5^PuRhGa5^UCoGa5 (forUCoGas the 
5/ states eventually cross above Ep)- 

At low energies, when > 0, all quasiparticle states 
are renormalized toward Ep, see the quasiparticle spectra in 
Figs. |3];c)-(f). In this energy region, is small, reflecting 
that quasiparticle states are coherent and itinerant. Above the 
peak in where < 0, quasiparticle states are pushed to 
higher energy. The lost spectral weight from the peak in is 
redistributed toward low energy near 1 eV in binding energy. 
A similar spectral weight redistribution occurs at the second 
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FIG. 3. (Color online) The computed momentum-averaged T>"{u)) 
and are plotted in (a) and (b), respectively All peak posi- 

tions in (y) in Fig. |2jb) are shifted to higher energy in T^" due to 
band-structure effects. Panels (c)-(f): The self-energy dressed angle- 
resolved spectral weight A(k,cj) is plotted along high-symmetry 
momentum directions. The peak-dip-hump feature is clearly evident 
in all spectra below 1 eV. 



peak (hump) in near 2 eV binding energy. As a result fur- 
ther pile-up of spectral weight occurs around 1 .0-1.5 eV, creat- 
ing new quasiparticle states due to electronic correlations. The 
quasiparticle states in this energy region are incoherent and 
fairly dispersionless, reflecting the dual aspect of the local- 
ized behavior of 5/ electrons. Qualitatively similar behavior 
was also found by using the LDA+DMFT method, however, 
with a weaker renormalization toward the Fermi level |20|. 

To compare our calculations with experiment, we compute 
the PES spectra as Ipes = (^(k, co'))^ni?(a;) (neglecting 
any matrix-element effects). We compare with available data 
for PuCoGas at 77 K ||9l shown by magenta diamonds in 
Fig. [4] Good quantitative agreement is evident. Near Ep ex- 
periment shows a broader feature than theory with less spec- 
tral weight, which may be related to experimental resolu- 
tion and theoretical approximations. The present calculation 
slightly underestimates the dip in the spectral weight, which 
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FIG. 4. (Color online) (a) Computed PES spectra for various com- 
pounds are compared with data for PuCoGas |9|. (b) Zoomed in 
view of (a) for PuCoGas spectrum. All theoretical spectra have been 
renormalied by the same scaling factor. 
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FIG. 5. (Color online) Experimental values of Tc and an impurity 
phase T* 1 3 1 are plotted as a function of the bare DOS at (theory) 
and compared with computed values of the spin-fluctuation coupling 
constant A and corresponding Sommerfeld coefficient 7. 



A from the energy derivative of S^ In the low-energy re- 
gion, we obtain S'(^k) ^ — -^^k = (1 ~ ^~^)^k- The cou- 
pling constant A follows the same material dependence as Tc 
across the series from PuCoIns^UCoGas with its maximum 
for PuCoGas . Although A is quite large for PuCoIns , its Tc 
is strongly suppressed probably due to competition with an 
impurity phase |3|. Our estimation of the fluctuation renor- 
malized Sommenfeld coefficient 7 follows Tc in Fig. [5] For 
PuCoGas, we find the renormalized 7 = 57 mJ/mol/K^, 
which is slightly less than the corresponding experimental 
value of 77 mJ/mol/K^ 1 1 1, suggesting room for phonon fluc- 
tuations of about Xep ~ 0.8, which is very close to the 
electron-phonon coupling constant deduced by first-principles 
calculations | 21 1. Note that our calculated coupling constant 
of A = 1.4 for PuCoGas is smaller than the calculated value 
of 2.5 obtained within the LDA+DMFT approximation ll22ll . 

In conclusion, we presented a first-principles based inter- 
mediate coupling model for calculating the multiband spin- 
fluctuation spectrum within the GW method. The presence of 
a strong spin-fluctuation peak in is found around 0.5 eV, 
which splits the electronic states into an itinerant coherent part 
close to and strongly localized incoherent states around 
1.0-1.5 eV. These results agree well with the experimental 
peak-dip-hump PES structure |9|. In addition, the isostruc- 
tural Pu-115 and UCoGas compounds (for UCoGas the 5/ 
electrons are moved above Ep) have qualitatively similar 
electronic band structure near Ep- This can be understood ap- 
proximately within a unified rigid-band filling scheme, which 
can account for band shifts through controlled hole doping. 
Finally, we calculated a spin-fluctuation coupling constant A 
of order unity. It follows the same materials dependence as T^ 
indicating that spin-fluctuation mediated pairing is a strong 
candidate for superconductivity in these materials. 
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stems from the neglect of orbital matrix-elements, charge and 
other fluctuations, as well as the quasiparticle approxima- 
tion in the self-consistency scheme of the calculation of the 
self-energy. The key result is that both the spectral weight 
loss at low energy and high energy are well captured by the 
spin-fluctuation model. As we move across the series from 
PuCoIns to UCoGas the spectral weight redistribution grad- 
ually decreases. This suggests that spin-fluctuations play a 
lesser role in UCoGas than in the isostructural Pu-115 com- 
pounds. 

Finally, we calculate the spin-fluctuation coupling constant 



SUPPLEMENTARY INFORMATION 

Intermediate coupling model— We calculate materials- 
specific first-principles electronic band structures, includ- 
ing spin-orbit coupling, within the framework of density 
functional theory in the generalized-gradient approximation 
(GGA) 1 12 |. We use the full-potential linearized augmented 
plane wave method of WIEN2k. The calculation is performed 
with 40 bands to capture the ±10 eV energy window of rele- 
vance around the Fermi level. The spectral representation of 
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the Green's function is constructed as 



Gspi}i,iuJn) ^ 



(3) 



Here is the eigenstate for /i^^ band projected on the i^^ 
orbital. The non-interacting susceptibiHty in the particle-hole 
channel represents joint density of states (JDOS), which can 
be calculated by convoluting the corresponding Green's func- 
tion over the entire Brillouin zone (BZ) to obtain (spin and 
charge bare susceptibility are the same in the paramagnetic 
ground state) ! 141 : 

2^ ^ 



k,n 



(4) 

Within the RPA, spin and charge channels become decoupled. 
(We ignore particle-particle as well as weaker charge fluctu- 
ation processes.) In the spin channel, the collective many- 
body corrections of the spin-fluctuation spectrum can be writ- 
ten in matrix representation: x = X^[i ~ ^sX^]~^- The in- 
teraction matrix Us is defined in the same basis consisting of 
intra-orbital inter-orbital V, Hund's couling J and pair- 
scattering J' terms |[T4HT6l . In the present calculation, we ne- 
glect the orbital overlap of eigenstates, i.e., we assume = 1 
when i = II. Such an approximation simplifies the calculation 
and becomes a diagonal matrix and J = J' = 0. 

Using the GW approximation, where G represents the 
Green's function and W is the interaction vertex, we write 
the spin-fluctuation interaction vertex following Ref . [il7i : 



^pqrs (Q? ^) — 



(5) 



The Feynmann-Dyson equation for the imaginary part of the 
self-energy in a multiband system with N sites is 

q rs "^-^ 

• ArsO^ + q, ^ - ^){nB{^) + 1 - nric^ - ^)) • (6) 



The quasiparticle spectral weight is defined by Ars (k, e) = 
— Im[Grs(k, e)] /tt. The quantities ub and np are the Bose- 
Einstein and Fermi-Dirac distribution functions, respectively, 
r is the vertex correction discussed later. For a more accu- 
rate calculation, one needs to account for the anisotropy in 
F(q, ft). In the present case, where the spin-fluctuation spec- 
trum is only weakly anisotropic (see Fig. 2), it is justified to 
replace the first term in Eq. ([6]) by a momentum- averaged 



spin-fluctuation function interaction, that is, (y{Q))^ = 
J d^qV{qj ft). This is equivalent to dropping the k 



(27r)3 

dependence of the self-energy, which greatly simplifies the 
numerical self-consistency loop. It then follows from Eq. ^ 
that at T = 0, the imaginary part of T>pq{k^ uj) reduces to 



KM 



dnT{Vpqrs{^)) Nrs{0J-ft), (7) 



for UJ > 0, where the density of states is given by A^rs(e) = 
Xlk "^rs (k, e) . (For u < 0, the only changes are that the upper 
limit of the integral is |a;| and the argument of Nrs is ft — \uj\, 
which is < 0.) We use Eq. ^ to compute the imaginary part 
of the self-energy from the first-principles band structure. The 
real part of the self-energy, E^g(co'), is obtained by using the 
Kramers-Kronig relationship. Finally, the self-energy dressed 
quasiparticle spectrum is determined by Dyson's equation: 



G-\k,uj) = G^\k,cu) - {t{k,uj))k. 



(8) 
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